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8 ! ABSTRACT 

Q 

We report on spectral and intensity variability analysis from a 
Chandra / ACIS-S observation of the transient, type-I X-ray bursting low- 
mass X-ray binary Cen X-4. The quiescent X-ray spectrum during this 
observation is statistically identical to one observed previously with BeppoSAX , 
and close, but not identical, to one observed previously with ASCA. The 
X-ray spectrum is best described as a pure Hydrogen atmosphere thermal 
spectrum plus a power-law component that dominates the spectrum above 2 
keV. The best-fit radius of the neutron star is r =12.9±2.6 (d/1.2 kpc) km if 
P H i the interstellar absorption is fixed at the value implied by the optical reddening. 

Allowing the interstellar absorption to be a free parameter yields r =19i 4 Q 

CO 



(d/1.2kpc) km (90% confidence). The thermal spectrum from the neutron 
star surface is inconsistent with a solar metallicity. We find a 3a upper-limit 
of root-mean-square variability < 18% (0.2-2.0 keV; 0.0001-1 Hz) during the 
^ ■ observation. On the other hand, the 0.5-10.0 keV luminosity decreased by 

40±8% in the 4.9 years between the ASCA and Chandra observations. This 
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variability can be attributed to the power-law component. Moreover, we limit 
the variation in thermal temperature to ^10% over these 4.9 years. The stability 
of the thermal temperature and emission area radius supports the interpretation 
that the quiescent thermal emission is due to the hot neutron star core. 

Subject headings: stars: neutron — stars: individual (Cen X— 4) 



Introduction 



Measuring both the mass (M) and the radius (R) of a neutron star (NS) would strongly 
constrain the nuclear equation of state. To distinguish between the competing models, 
several measurements of a few percent accuracy in M and R are required. While the masses 
of several NSs have been constrained to 10% or better (see [l'horsett fc Chakrabarty 1999 



for a recent review) from pulsar timing of Doppler shifts, measuring the NS radii has 
proven to be difficult. There are at least five ways to measure R using X-ray emission from 
accreting NSs: (1) the spectral evolution of radius-expansion type-I X-ray bursts; (2) the 
measurements of M/R from the gravitational red-shift of metallic spectral lines during 
type I X-ray bursts or 7-ray lines during accretion; (3) inferring constraints from NS kHz 
quasi-periodic oscillations; (4) pressure broadening and red-shift of photospheric metal 
lines; and (5) spectral analysis of transiently accreting NSs in quiescence. 

Efforts to measure the NS radius through the spectral evolution of radius-expansion 
type-I X-ray bursts (for a review, see [Lewin et al. 1993| ) based on theoretical non-Planckian 
spectra ( [Van Paradijs 1982 ; London et al. 1986 ; Pavlov et al. 1991 ; Made] 1991 ; 



Titarchuk 1994|) have been somewhat successful. Performing this measurement with 
the emission from a type-I X-ray burst has the advantage that the observed luminosity 
originates from the NS photosphere and not the surrounding accretion disk. The measured 
NS radii range from 6 to 15 km. The reliability of these measurements is limited by 
systematic uncertainties in the emergent spectrum, the fraction of the NS surface area 
involved in the burst and the elemental composition of the photosphere, as well as the 
distance to the NS ( |Van Paradijs fc Lewin 1987| ; pamen et al. 1989| ; |Damen et al. 1990| ) . 



Attempts to measure NS photospheric metal lines as expected from model atmospheres 
( poster et al. 1987|) have had mixed success. Absorption lines observed in the tails of type-I 
X-ray bursts ([Waki et al. 19841 ; [Nakamura et al. 1988| ; |Magnier et al. 1989| ) have not been 
confirmed through repeated observation with more sensitive instrumentation, resulting in 
metallicity limits of Z < 0AZ Q ( pay et al. 1992| ). Moreover, the interpretation of these lines 
as due to absorption in the photosphere is in doubt ( |Madej 1990| ). Repeated observations 
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of such lines, and their identification with a known transition would measure M/R. 
Gamma-ray lines from the accretion process are also a possibility flBildsten et al. 1992| ; 
Bildsten et al. 1993| ), though the levels of emission are still below the ability of current 
instruments. 

Observations of kilo-Hz quasi-periodic oscillations from accreting neutron stars, 
interpreted as the orbital period at or above the marginally stable orbit ( |Kluznial 
|et al. 1990| ; [Kaaret et al. 1997] [Kluzniak 1998| ) have been used to constrain both M and R 
using observations from eight accreting neutron stars, finding M ^2M (|Zhang et al. 1997| ). 
The NS radius is then less than the marginally stable orbit for a mass of this size (<18 km 
for a non- rotating star). 

Another method of measuring R for a slowly rotating NS is the approach described 
by Paerels ( |1997|) , with metal lines from an emergent NS X-ray spectrum. Measuring the 
pressure broadening (oc M/R 2 ) and photospheric red-shift (oc M/R) yields both M and R. 
This method is intriguing because it is independent of distance, which can be uncertain by 
up to a factor of 2 in the isolated transient NSs, although it can be determined to ~5% 
for the low-luminosity X-ray sources in globular clusters ( [Hertz fc Grindlay 1983j ; |Verbunt 



et al. 198~3 ). Once several metal lines are identified in the photospheric spectrum, this 
method is very likely to provide a number of M and R measurements. 

To these approaches has recently been added a new means of measuring R with X-ray 
spectroscopic observations of transiently accreting, low magnetic-field (<10 10 G) NSs in 
quiescence (|Brown et al. 1998| |Rut ledge et al. 1999] ; [Rutledge et al. 2000]) . Until this 
work, the spectral observations had been interpreted as black-body radiation ( |Van Paradijs 
et al. 1987| ; |Garcia 199^ ; |Verbunt et al. 1995| ; |Asai et al. UMb] ; |Asai et al. I^B] |Uampana| 
et al. 1998bj ; |Garcia fc Gallanan 1999| ; |Uampana et al. 2000| ), yielding effective temperatures 
of /cToo=0.2-0.3 keV. For most sources, the thermal component has a black-body emission 
area radius much smaller than a NS (<;1 km). 



As described by Brown, Bildsten & Rutledge ( |1998j , BBR98 hereafter), a black- 
body spectrum is not appropriate for weakly magnetic 7 transiently accreting NSs with 
/cT e fj < 5xl0 6 i^. At accretion rates <C 10 12 g s -1 the accreting metals gravitationally 
settle faster than they are supplied, and the atmosphere is nearly pure hydrogen ( |Bildsten 



et al. 19§2] ). The dominant opacity is free- free absorption (oc v ) which results in a 
spectrum in which higher energy photons escape from deeper in the NS atmosphere, where 
the temperature is higher (Pavlov fc Shibanov 1978] ; |Rajagopal fc Romani 19961 ; |Zavlin 



Jn low magnetic field NSs (<10 10 G), the magnetic field plays no role in determining the opacity of the 



atmosphere at energies above 0.1 keV, and thus can be neglected ( Zavlin et al. 1996 ) 
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et al. 1996). This hardens the spectrum which, (mis-) interpreted as a black-body, results 
in systematically higher temperatures and lower emission area radii. When the ASCA and 
ROSAT observations are fit with the more appropriate thermal H atmosphere spectrum, 
the emission area radii are ~10 km flRutlcdgc et al. 1999 ; Rutlcdge et al. 2000 ), confirming 
that this emission is mostly thermal. These objects have thus become the focus of our 
efforts to measure the NS surface area. 

If the emission area radius of the thermal part of the quiescent X-ray spectrum is the 
NS radius, it should be constant from observation to observation. In addition, if accretion 
onto the NS makes no contribution to the quiescent luminosity, the temperature should 
also be constant on timescales -ClO 6 yrs, the core cooling timescale (see |Ushomirsky & 
Rutledge 2001| for exceptions). In the most sensitive observations, an additional power-law 
component is observed above a few keV ( [Asai et al. 1996b] ; |Asai et al. 1998| ; |Campana 
et al. 1998b| ; Campana et al. 2000 ). If produced through interaction between the accretion 
and magnetosphere QCampana et al. 1998b| ), it might vary on short (~msec) timescales, 
and perhaps show a pulsation at the NS spin frequency. 

In this paper, we present X-ray spectroscopic analysis of Cen X— 4 in quiescence. 
Cen X— 4 is a transient, type-I X-ray bursting low- mass X-ray binary (for reviews on 
transients, see |Tanaka fc Lewin 1995| ; |Tanaka fc Shibazaki 19961 ; |Chen et al. 19971 ). Its 
distance is estimated to be 1.2±0.3 kpc on the basis of two observed radius-expansion 
bursts flChevalier et al. 1989|) ; formally, this is a distance upper-limit. It is in a P or b = 0.629 
day binary with a K dwarf QChevalier et al. 1989| ; [McClintock fc Remillard 1990|) , and has 
been observed in outburst twice (1969 and 1979; see discussion in § |3J]). We compare the 



observed X-ray spectrum with a number of models, in particular with the H atmosphere 
model. We also compare the spectrum with spectra of the same source in archived 
observations. We place limits on changes in the thermal spectrum of the NS on year long 
timescales. In §Q, we present the analysis of the Chandra data. In §2.3] , we compare the 
Chandra data with past observations, and characterize the spectral differences between this 
and past X-ray observations. We discuss these results and their implications for this and 
other sources and conclude in SBL 



2. Observation and Analysis 

The observation occurred 23 Jun 2000 01:23:01-04:32:54 TT for a total exposure time 
of 9561.4 seconds with Chandra / ACIS-S3 (backside illuminated). The source position was 
aimed 4' off-axis, with a 1/8 subarray used with 0.4 sec exposures. The time between 
successive frames was 0.44104 sec, which gives ~ 10% dead-time. We analyzed data from 
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the pre-processed Level 2 FITS data provided with the standard data products. The X-ray 
source appears at the known optical position of Cen X-4 (±1 "; |Canizares et al. 1980 



). Only one other point source (previously unknown) was found in the field with the 
Chandra Interactive Analysis of Observations Software (CIAO) point-source detection tool 
celldetectf\> with a S/N of 3.5, roughly 1' away. 

We extracted data from a circle 10 pixels in radius about the Cen X— 4 X-ray source 
position, and background from an annulus with inner and outer radius of 13 and 50 pixels, 
respectively. There were a total of 2714 good counts in the source region, and 575 in the 
background region; we expect 25 background counts in the source region (~1% of the total 
counts in the source region). 



2.1. Chandra Spectral Analysis 

The temperature of the focal-plane instruments has been decreased during Chandra's 
lifetime; this alters the energy response of the ACIS-S chips. The data were taken with a 
focal-plane chip temperature of — 120C, and we used the corresponding response files for 
this temperature, according to the standard ACIS-S analysis. 

We binned the photons between 0.5 and 1.5 keV into bins of width ^130 eV wide, 
comparable to the energy resolution (FWHM) in the BI-S3 detector. Above 1.5 keV, we 
binned the data so that there were 40 counts per bin (wider than the spectral resolution). 
We used spectral data in the energy range 0.5-10 keV. While there are significant counts 
below 0.5 keV, the ACIS-S energy calibration is presently not reliable below this energy; 
we will re-examine this analysis when the energy calibration is refined below 0.5 keV. In 
addition, the current Chandra response files for the ACIS-S-BI underestimate the detector 
area by up to 20% near 0.5 keV (N. Schulz, priv. comm.). To account for the calibration 
uncertainties, we included a 25% systematic uncertainty in the 0.5-0.6 keV energy range, 
and a 5% systematic uncertainty in 0.6-0.7 keV energy range. This is in addition to a 4% 
systematic uncertainty across all energies that accounts for other calibration uncertainties. 



We fit the data in XSPEC vll ( [Arnaud 1996[ ), with several spectral models (powerlaw 



H atmosphere, blackbody, Raymond-Smith, or thermal bremsstrahlung) , all with galactic 
absorption (Ah) as a free parameter, except where noted. The H atmosphere spectrum is 
not a standard XSPEC model, but has been described and calculated elsewhere ( [Zavlin 



et al. 1996|) . The hydrogen column density for Cen X— 4 was estimated during quiescence 



8 available at 



http:/ /asc. harvard.edu 
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to be iVH,22<0.20 ( |Asai et al. 1996a| ; A^ H =10 22 iVH i 22 cm 2 ); its optical reddening is 



E(B — V) = 0.1 ( |Blair et al. 1984j ). These values are consistent with an approximate optical 
reddening/hydrogen column density ratio found from observing the halos of X-ray sources 
( |Gorenstein 1975 ; Predehl fc Schmitt 1995 ), which imply an equivalent hydrogen column 



density of Ah,22 =0.055. 

No single component model fit the data acceptably (prob <C 10 -6 ). In the black-body 
and H atmosphere fits, the high-energy (>2 keV) powerlaw spectral component reported 
previously ( |Asai et al. 1996b] ; |Campana et al. 2000|) is apparent. In general, the single 



component models failed due to the presence of this high-energy component. We then fit 
the data with two-component models, where the second component is a power-law that 
accounts for the >2 keV emission. The best-fit absorbed H atmosphere spectrum has 
parameters listed in Table [I]. We also provide the best fit parameters with Ah held fixed at 
the optically implied value (Ah, 22=0. 055), and the resulting intrinsic (that is, unabsorbed) 
spectrum in Figure [p. 

While a double power-law model was statistically acceptable, the low energy component 
has a steep power-law slope (photon index a=6.2tj;g, 90% confidence; all uncertainties 
and upper-limits are 90% confidence unless otherwise stated) and considerably higher 
absorption than observed previously from this source or implied by optical observations 
(^h,22=0.501q; 2 q). We reject the model on this basis. 



A Raymond-Smith plasma model ( [Raymond fe Smith 1977 ) with an underlying 



power law is also statistically acceptable, however, with an abundance limited to 
Z <1.4x 10~ 3 Z so i. In addition, it has been shown previously that the X-ray-to-optical 
flux ratio is substantially greater than is typical for stellar coronal sources flBildsten fe| 
Rutledge 2000|). We therefore reject a model which ascribes this emission to the corona of the 



stellar companion. Though not physically motivated, a thermal bremsstrahlung spectrum 
acceptably fits the data, with a volume emission measure of / n e nidV = 2.3xl0 56 (ci/1.2kpc) 
cm -3 and a temperature of 0.34 keV. 

To enable comparison with previous work, we also fit the data with an absorbed 
blackbody and powerlaw. The best fit with a fixed A^h,22=0.26 QCampana et al. 2000| ) is 



9 While the metallicity of the companion has not been observationally constrained, Martin et al.(1994), 
while investigating Li abundance, found that the Ca I A6717 line is consistent with solar metallicity, and 
that an assumed solar Fe abundance in their spectral fits to the Li I doublet is consistent with the data 



In optical spectral studies of the companion ( Cowley et al. 1988 ; |Chevalier et al. 1989 ; McClintock fe 



Rcmillard 1990 ; McClintock fe Remillard 2000 ), while metal absorption features are observed, the metallicity 



of the companion is not estimated; however, it does not appear to be so substantially sub-solar. 
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rejected (x^=3.1/10 dof; prob=5xl0 4 ). With all parameters free, we find A^h,22<0.071 
(90%); fcT eff =0.175t£^; r=l.l±g;| (d/1.2 kpc) km; a =1.2+°;| ( x 2=0.75/9 dof; prob=0.67). 



2.2. Solar Metallicity Atmosphere 

We fit the Chandra data with a solar metallicity atmospheric model for a 10 km, 1.4M Q 
NS, with relative abundances given by Grevesse & Noels ( |1993| ). A galactic absorption and 
solar metallicity atmosphere model failed to fit the spectrum (x1= 18/11 dof), due to the 
high energy (>2 keV) excess attributed above to the power-law component, and to spectral 
curvature below 1 keV. With an additional power-law component, the best-fit spectrum 
is still unacceptable (prob = 3xl0 -8 ), largely due to the Fe L edge near 0.534 keV (for a 
1.4M Q , 10 km NS). While the model atmospheres assume no rotation and are static, the 
major discrepancy between the model and data is this absorption edge, which will not be 
significantly altered by these effects. 

A solar metallicity NS atmosphere spectrum is therefore strongly rejected. More 
detailed investigation regarding sub-solar metallicity spectra is in progress, and will be 
presented in forthcoming work. It would be highly useful to have response matrices refined 
down to 0.3 keV for this work. These will be valuable limits, as accretion at a rate high 
enough to explain the quiescent thermal emission would enrich the atmospheric metal 
content to a detectable level (BBR98). 



2.3. Comparison between the Chandra, BeppoSAX and ASCA X-ray spectra 

Cen X— 4 has been detected on 6 occasions in quiescence, with Einstein and EXOSAT 
( |Van Paradijs et al. 1987|) , twice with ASCA ( |Asai et al. 1996b] ; [Asai et al. 1998Q , once with 



ROSAT /HRI (|Campana et al. 1997| ) and once with BeppoSAX (|Campana et al. 2000| ). The 
measured luminosities and blackbody spectral parameters from the published observations 
are listed in Table 0. 

In this section, we compare the Chandra X-ray spectrum with that previously obtained 
with BeppoSAX flCampana et al. 2000|) and with the ^^CA/SIS+GIS spectrum ([Asai 



et al. 1996a| ; |Asai et al. 19981 ; [Rutledge et al. 1999Q . We do not compare with one other 



existing ASCA/G1S observation as it was performed largely with the GIS, with lower S/N 
than the first ASCA observation, and is comparable to the first ASCA observation ( |Asai| 
et al. 19981) . We have previously re-analyzed the ASCA observation ( [Rutledge et al. 1999|) 



and use the same resulting spectrum here. 
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2.3.1. Re- Analysis 0/ BeppoSAX observation 

Campana et al.(2000; COO hereafter) interpreted the quiescent X-ray spectrum 
(measured with ASCA and BeppoSAX) of Cen X— 4 as a black-body, and found a mean 
blackbody radius 3.1 km for the thermal component (the best-fit blackbody radius with 
BeppoSAX data alone was r=10li° 8 8 km )- However, N n was held fixed at the best-fit value 
of -Wh,22=0.26, which is higher than the optically implied value. We undertook a reanalysis 
of the BeppoSAX observation of Cen X— 4, previously analyzed by COO. There are two 
major differences between our spectral analysis of the LECS data (0.3-2 keV) and theirs: 



1. For the LECS background data, COO used the standard "blank fields" background 
which are taken from several "blank" areas of the sky (Parmar et al. 1999|) . However, 



Cen X— 4 is located in a region of the sky with higher than average low-energy 
background (/ = 332.24 deg, b = +23.88 deg), according to the ROSAT/All-Sky- 
Survey Soft X-ray Background maps (~0.75 keV; [Snowden et al. 19"97| , Fig 6b). Thus, 
the "blank fields" background underestimates soft X-ray background in the Cen X— 4 
observation. We used an alternative method of background estimation described by 
Parmar et al. ( the "annulus" method; |1999[), in which the background intensity is 



assumed to be proportional to the countrate in an annulus in the LECS detector, 
with a spectrum identical to that found with the "blank fields" method. By this 
method, the 0.3-1.0 keV background counts account for 60% of the total counts in 
the Cen X— 4 source region, whereas they only account for 34% in the "blank fields" 
method; the counting uncertainty in the source region is 10%, and therefore this 
difference in background countrate is significant. As the "annulus" method scales 
with the local background countrate, it is a more realistic estimation of the local 
low-energy (<2 keV) background, and so we adopt it. 

2. For the LECS data, we used a larger extraction radius than COO (8 ' vs. 4'); the 
LECS point-spread-function becomes broad toward the lower energies (80% encircled 
energy at 1.5 keV is 3'.0, and at 0.28 keV is 6'.1; the 95% radius is 5'. 5 and 8'. 5 at 
1.5 and 0.28 keV, respectively). The smaller extraction radius used by COO may have 
biased the resulting Ah upwards. This collected 647 counts from the LECS in 21524 
sec (much greater than the 233 counts found by COO). In the 0.3-1.0 keV, there are 84 
counts in the source region. We estimate that using an 8' radius instead of a 4' radius 
increases the observed source countrate by 30%, which is larger than the Poisson 
uncertainty of 10%. 



For the MECS data, we used a 4' extraction radius, which extracted 651 counts 
(comparable to the 632 found by COO). We used the standard MECS background file 
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(MECS23_bkg.evt) extracting background counts using the same region as for the source 
fields. The source background subtracted countrates are (5.3±1.6) x 10~ 3 c/s (0.1-3.1 keV) 
and (3.3±0.4)xl(T 3 c/s (1.7-9.0 keV) in the LECS and MECS instruments, respectively. 
The LECS data has lower S/N than reported by COO due to the higher background; we find 
a slightly higher MECS countrate than found previously. 

For comparison with previous work, we fit the BeppoSAX LECS+MECS data 
with an absorbed black-body plus power-law spectrum. We obtained an acceptable fit 
(x^=0.65/l dof; prob=0.42), with the following best-fit parameters: A" H ,22<0.95, a=1.6io;|, 
kT eS =0. 10io!o4 keV. We were unable to put any reasonable limits on the BB radius (>0.4 
km, unbounded from above), as the best-fit spectrum merely would increase the column 
density to compensate for larger and larger areas. When we hold the column density fixed 
at the best-fit value, we find rBB=7-l% 3 km; when we hold the column density fixed at the 
optically implied value (A"h,22=0.055), we find r BB =2lll s km. These values are comparable 
to those found with ASCA (|Asai et al. 1996b| ; |Rutledge et al. 1999|) , and below those 



found by COO. We attribute this difference to our more accurate background subtraction 
and larger extraction radius for the LECS data. The unabsorbed total flux is 2.7xl0~ 12 
erg cm -2 s -1 (4.7xl0 32 ergs -1 , 0.5-10.0 keV); the unabsorbed flux of the BB component is 
2.2xl0~ 12 ergcm~ 2 s _1 (0.5-10.0 keV), and of the power-law component it is 0.44xl0~ 12 
ergcm _2 s _1 (0.5-10.0 keV). As we find in the next section, the BeppoSAX and Chandra 
spectra (and their corresponding fluxes) are, within statistics, the same. 



2.3.2. Joint Spectral Fitting of Chandra, BeppoSAX and ASCA X-ray Spectra 

We performed a joint spectral fit of the Chandra, ASCA, and BeppoSAX data. We 
assumed a 4% systematic uncertainty in the detector responses for all instruments, and 
held N H fixed at its optically implied value. We used an absorbed H atmosphere+powerlaw 
spectrum. The best spectral fit rejects a single model to account for all three observed 
spectra (x 2 ,=3.3/53 dof; prob=5x 10~ 15 ). A joint Chandra/ BeppoSAX fit provided a 
statistically acceptable fit, indicating that the Chandra and BeppoSAX spectra are 
statistically identical. A joint Chandra/ ASCA fit was not acceptable. The ASCA data are, 
by themselves, acceptably fit by this same assumed spectrum ( |Rutledge et al. 1999j ). 

We jointly fit the Chandra and ASCA spectra, permitting one of the five parameters 
(kT, r, a, and power-law normalization N a ) to be different between the two spectra, each 
in turn. None of these fits are formally statistically acceptable. However, while a change 
in a is soundly rejected (prob<4x 10~ 10 ), and a change in r is of low probability (prob 
=0.001), changes in the N a or kT are found to be marginal (both prob = 0.01). This is in 
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the range where systematic uncertainty in the detector responses (in particular, the area 
as a function of energy) become important. For example, if we increase the systematic 
uncertainty from 4% to 8% (which, for example, could be due to a systematic offset in the 
absolute flux calibrations between the two detectors), the best fits for a changing N a or kT 
become statistically acceptable. 

Thus, while the Chandra and ASCA spectra are significantly different, we cannot 
unequivocally state whether this difference is due to a change in the N a or the kT, a 
combination of these with other parameters, or a systematic difference between the absolute 
calibrations of the two instruments. However, we provide the best-fit spectral parameters 
in Table |3|, for a changing power-law normalization and thermal temperature. These then 
serve as an upper-limit for variability in these parameters between the two observations. 
Between the Chandra and ASCA epochs, then, we find that the thermal temperature 
decreased (at most) from kT =0.085 to 0.077 keV; or the power-law normalization decreased 
from 15.2xl0~ 5 to 5.9xl0 -5 phot cm -2 s" 1 at 1 keV. Moreover, the interpretation of these 
spectra impacts the best-fit power-law slope, which is a=1.0 if it is the temperature which 
changes, or a— 1.7 if it is the normalization which changes. 

Finally, we fit a spectral model to the Chandra, BeppoSAX , and ASCA data assuming 
that the power-law slope and flux changes, but that the thermal component does not. 
The resulting best-fit parameters and uncertainties are shown in Table The best-fit 
is statistically acceptable. The power-law component (taken to the be the same for the 
Chandra and BeppoSAX observations) is flatter in the Chandra+ BeppoSAX observations 
than during the ASCA observation (1.0±0.4 vs. 1.9±0.3), and has a lower flux, by about a 
factor of two. 

We also fit these spectra assuming that the thermal r and kT vary while the power-law 
component does not (although this is not theoretically motivated). An acceptable fit is 
found (x^=l-33/52 dof; prob=0.06), with Chandm+BeppoSAX values of (r=13.4jr 2 ;4 km , 
kT=0.07At° f t keV) and ASCA values of (r=7.1±H km, fcT=0.10lffikeV). 



2.4. Intensity Variability 

Variability in the luminosity of Cen X— 4 has been observed over timescales from 
days to years. Van Paradijs et al. (|1987|, JVP87 hereafter) detected Cen X— 4 using 
EXOSAT/LE1 with the CMA instrument and LEX AN 3000 A filter flTaylor et al. 1981 



with a countrate of (5.4±1.2)xl0 3 c/s. JVP87 assumed a thermal bremsstrahlung 
spectrum with temperature between 1 and 5 keV and a column density iYn 22 =0.066, which, 
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for this countrate, corresponded to a flux at earth of (2.2-5.6) xlO -12 erg s -1 cm" 2 . Using 
the same assumptions with an earlier Einstein/lPC observation, they found a flux at earth 
between (1.1-1.5) xlO -12 erg s" 1 cm -2 , concluding that the X-ray luminosity must have 
increased between the Einstein and EXOSAT observations by 50-500%. Campana et al. 
( |1997|) measured a flux variation of a factor of 3 in <4 days using ROSAT /HRI. 

We started by looking for variability during our observation. We first separated the 
data into two energy bins: pulse-invariant (PI) bins 2-136 (0.2-2 keV) and PI bins >136 
(>2 keV). We included PI bins in which the area and energy calibration is presently not 
reliable; while these are not useful (presently) for flux and spectral calculations, they can 
be used for investigations in variability. 

There were a total of N phot = 2553 counts in PI bins 2-136 («0.2-2 keV). We used 23903 
time bins which were 0.44104 sec in length which each had exposure of only 0.40 sec, for a 
total observation time of T Q b s = 10542 sec. We performed a Fourier transform, producing a 
power density spectrum (PDS) of 11951 frequency bins across the 0.0001-1.13 Hz frequency 
range, normalized according to Leahy et al. ( |1983| ). We logarithmically rebinned this data, 
and fit the PDS with functions -P(z^), to extract the integrated r.m.s power. The PDS is 
acceptably described as counting noise (constant power of 2.0; x^=l-41/16 dof; prob=0.13); 
we therefore observe no variability in this data. Using a power-law distribution with the 
slope of a = 1 held fixed, with an underlying Poisson level (P = 2) also held fixed, we find 
a 3cr upper limit on the root-mean-square variability of <18% (0.0001-1 Hz). For a flat 
power-law (a = 0), the 3a upper limit is <10%. 

There are a total of 151 counts in PI bins >136 (^2.0 keV) in this observation. We 
performed an identical PDS analysis as for the low energy counts, although in the fit to the 
PDS we held the exponent fixed at a=1.0 (producing an acceptable fit; xt~ 1-46/15 dof). 
The 3d upper-limit to the RMS variability >2.0 keV is <50% (0.0001-1 Hz). 

We also looked for longer term variability by investigating the hypothesis that the 
Chandra and ASCA spectra were identical in all parameters, but different in absolute 
normalization. We found that the spectra were describable in this way (% 2 = 1.05/47 dof; 
prob=0.37), with the ASCA spectrum a factor 1.65±0.12 (90% uncertainty) more luminous 
than the Chandra spectrum. We interpret this as fitting a well-constrained spectrum from 
Chandra (signal-to-noise ratio S/N^100) to the low S/N data of ASCA (S/N«4). This 
provides a rough measure of the luminosity difference between the ASCA and Chandra 
observations - a decrease of 40±8% over 4.9 yr. 
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3. Summary, Conclusions and Implications 

We have analyzed the X-ray spectrum and intensity variability of Cen X— 4 observed 
with Chandra, and compared it with earlier observations by ASCA and BeppoSAX . The 
Chandra X-ray spectrum is inconsistent with all single component models we applied, 
requiring at least a two component model, with a power-law component that dominates the 
spectrum above 2 keV and a softer (most likely thermal) component which dominates the 
spectrum below this energy. The intensity variability during the Chandra observation is 
<18% (3a) on timescales between 1-10,000 sec in the 0.2-2.0 keV energy band dominated 
by the thermal component. The limit on variability is weaker above 2 keV (<50% r.m.s., 
3d). 

Cen X-4's luminosity decreased between the ASCA and Chandra observations (4.9 
yrs) by 40±8% (0.5-10.0 keV). We attribute this variation to changes in the power-law 
component. This is comparable to the factor of ~3 decrease found over a few days by 
Campana et al. (|1997|) . We cannot statistically exclude that the power-law component 
remains constant, and the thermal component varies. Thus, Cen X— 4 is variable in 
quiescence on timescales longer than hours. 



3.1. The Thermal Component 

Rutledge et al. ( |1999| ; |2000| ) have shown that the soft emission seen in all NS transients 
in quiescence is best explained as thermal emission from a pure hydrogen NS atmosphere. 
We have found the same in this Chandra observation. Indeed, comparing X-ray spectra 
taken with Chandra, ASCA and BeppoSAX at different times over 4.9 years, we find 
that this thermal component is consistent with being constant, with a best-fit radius of 
r=12.9±2.6 (<i/1.2kpc) km, and a temperature fcT=0.076±0.007 keV. These results improve 
the precision of the quiescent spectrum by a factor of ~3 over previous work flRutlcdgc 



et al. 19991) . We place an upper-limit on the amount of temperature variability across 4.9 



years of ^10% (Table |3|). The agreement between this radius and that measured during 
Type I X-ray bursts from other systems is strong confirmation that the emitting area is the 
NS surface. 

Two energy sources for the thermal emission have been discussed: accretion in 
quiescence at a low-rate ([Van Paradijs et al. 1987] ; |Menou et al. 1999| ) and re-emission of 



heat deposited into the crust during the large accretion events (BBR98). The required 
accretion rate for the bolometric luminosity of 7 x 10 32 erg s _1 is M q « 6 x 10~ 14 M Q yr _1 , 
adequate to keep the metal content in the atmosphere comparable or larger than that in 
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the accreting material ( |Bildsten et al. 1992 ). If the metal content in the accreting material 
is very sub-solar (that is, effectively pure hydrogen), accretion at these rates would give 
thermal emission much like we observe ( |Zampieri et al. 1995[ ). However, solar metallicity 
accretion is likely ruled out, largely due to the lack of absorption from the Fe L edge near 
0.5 keV. Tighter constraints on the photospheric metallicity will be the subject of future 
work with these data and will allow us to more thoroughly constrain the active accretion 
hypothesis, which does not have a specific way of predicting M q . 

The other possible mechanism is re-emission of heat deposited in the inner crust due to 
pycnonuclear reactions, electron captures and neutron emissions ( |Haensel fe Zdunik 1990| ) 
during accretion events (BBR98). The layers where this heat is deposited are in close 
thermal contact with the NS core. The reactions then heat the core to a temperature of ~10 8 
K over ~10 4 yr (BBR98; |Colpi et al. 2001] ). The NS reaches an equilibrium where the time 
averaged nuclear heating equals the quiescent thermal emission, L q « (1 MeV/nucleon)(M), 
implying a time averaged accretion rate for this NS of (M) « 1O _11 M yr -1 . 

We can compare this (M) to that implied by the outbursts. The outburst in 1969 
had a total fluence of ~ 3 ergs cm" 2 ( |Chen et al. 1997] ). There has not been another large 
outburst of this magnitude recorded since; the small outburst in 1979 had a fluence two 
orders of magnitude smaller. It is uncertain to measure the time-averaged accretion rate 
on the basis of a single outburst. If a large outburst like that in 1969 occurs every 100 
years, the time-averaged accretion flux would be ~ 10~ 9 ergs cm" 2 s -1 , which at the 1.2 
kpc distance gives (M) ~ lO~ 11 M yr~ 1 , comparable to that implied by the level of the 
quiescent thermal emission. However, this is a low accretion rate for a binary at this long 
orbital period, where the likely driver of mass transfer is nuclear evolution of the binary 
(see |Webbink et al. 19"83|) . So, this is an outstanding puzzle for this source. If outbursts like 
that seen in 1969 occur more frequently than 100 years, the time averaged accretion rate 
would be higher and in conflict with that implied by the quiescent thermal emission through 
the mechanism of BBR98, unless the fraction of deposited energy re-emitted as photons is 
much less than unity. Possible solutions to this, involving enhanced neutrino emission from 
the core, have been proposed by Colpi et al. (|2001|) and Ushomirsky & Rutledge (|2001|) . 

If the thermal emission is due to a hot neutron star core, then the measured effective 
temperature of 0.076 ± 0.007 keV tells us the internal NS temperature. For example, if 
the outer layer consists of light elements to a density 6 x 10 5 gcm -3 , then using the fit 
of Potekhin et al. ( |1997[ ), we find T core = (5.ll];; 2 )xl0 7 K. For these core temperatures, 
the modified URCA neutrino luminosity is orders of magnitude less than the photon 
luminosity. The NS core would not be in a thermal steady-state if pion condensation 
occurs, as the resultant neutrino luminosity QUmeda et al. 19941) would be much larger than 
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the heating supplied by reactions in the deep crust. As a result, if some enhanced cooling 
mechanism were operating in the core, then we would expect T core ^ 10 7 K and the crust and 
surface temperature should be decreasing on a timescale of ~ 1 yr (BBR98; |Ushomirsky fc 
Rutledge 200l| ). 



The lack of luminosity variability in general, and the stability of the thermal component 
over a ~5 year timescale supports the model proposed by BBR98 for the thermal component. 
It remains to be unequivocally demonstrated, however, that the variability observed is 
restricted to the power-law component, and that the variability of the temperature of 
the thermal component is limited to ^1%, as expected from variations in the T core from 
outburst to outburst flColpi et al. 200 1[ ). 



3.2. The Power-Law Spectral Component 

The origin of the hard power-law component remains unclear. Many ideas have been 
put forward, ranging from emission due to an active pulsar wind colliding with the accretion 
disk, x-ray emission from a turned-on radio pulsar ( [Stella et al. 1994] ) and accretion onto 



the magnetosphere (for a review, see Campana et al. 1998a| ) . 



There is no direct knowledge of either the NS spin or magnetic field in Cen X-4, though 
the presence of Type I bursts points to B < 10 10 G. None of the above models make specific 
predictions about either the level of emission or its spectral shape which explain the present 
observations. 

Moreover, the limitation to <1% rms of coherent intensity pulsations in the transiently 
accreting NS Aql X-l (another transient, type-I X-ray bursting source) following a rapid 
decrease in flux ( |Chandler fc Rutledge 2000|) which was interpreted as due to magnetic 
inhibition of the accretion flow ( |Campana et al. 1998b| ; ^hang et al. 1998| ) does not support 
the idea that accretion onto the compact object is ultimately restricted by a "propeller 
effect" , favoring instead an interpretation of the end of transient outbursts as due to an 
end of the disk instability, such as occurs in dwarf novae ( [King fc Ritter 1998| ). In fact, 



abrupt declines in dwarf novae outburst fluxes are observed, and are explained within the 
disk instability model flCannizzo 1994 ), without invoking magnetic inhibition. The absence, 



then, of a magnetosphere in transiently accreting NSs would preclude attributing the 
power-law to magnetospheric accretion, and energy production would have to be attributed 
to another site within the NS binary system. 

Spectral comparison between the ASCA and Chandra observation indicate that a 
change in (only) the power-law normalization is marginally acceptable, and therefore that a 
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change in the spectral slope (from a = 1.9±0.3 to a = 1.0±0.4) is also required. It is unclear 
how a change in the spectral slope could be produced while limiting the luminosity decrease 
to less than a factor of three in any of the proposed emission mechanisms ( |Campana 
et al. I998a| ). 



There might be some help in understanding the power-law component from other 
wavelength bands, where it was found recently that, in quiescence, the energy density (uF v 
is nearly flat from UV through the X-ray energy range, in marked contrast to A0620-00, 
in which the energy density falls in the UV band ( [McClintock fc Remillard 2000|) . This 



flat-spectrum has been interpreted as due to a shock at the splash-point of accretion 
from the companion flMenou fc McClintock 2001[ ). However, while we also find a flat 



spectrum for the ASCA observation (photon power-law slope of a = 1.9±0.4) we find an 
increasing energy density during the Chandra observation (a = 1.0±0.4). Clearly, further 
observational study of the power-law component is required to understand even its most 
basic aspects. 



3.3. Future Work on Neutron Star Radii from Thermal Emission 

The dominant systematic uncertainty in the measurement of the radius of the neutron 
star in Cen X— 4 is the distance to the object. This will be dramatically improved by the 
Space Interferometric Mission, which can measure the parallax of the 18. m 5 counterpart to 
~ 10 yuasec, or <1% distance uncertainty at 1.2 kpc. This would make the dominant 
uncertainty in the measured thermal emission radius, followed by any systematic uncertainty 
in the modelled spectrum. For example, a planned 50 ksec XMM observation will permit 
constraint of r/D to ~2% while simultaneously measuring N^. 

It has been suggested that some fraction of the low-luminosity X-ray sources in 
globular clusters are transiently accreting NSs in quiescence ( [Verbunt et al. 1984] ; |Rut ledge" 



|et al. 2000| ) . These make excellent targets for NS radius measurements, as there are 
multiple objects per observing field, all at the same distance and interstellar column density 
(BBR98). Chandra imaging observations of globular clusters ( prindlay et al. 2000| ; |Pooley| 
|et al. 2 0"0(]| ) indicate source densities in excess of several per square arcmin (down to 
luminosities of Lx ~10 30 ergs -1 ) with over 100 X-ray sources in the core of 47 Tuc alone. 
It is still unanswered as to what fraction of these are quiescent NS's. The high source 
density will complicate X-ray spectroscopy with XMM or Con-X (with angular resolution 
~15"), making Chandra imaging spectroscopy the best way to pursue this science. 
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Fig. 1. — The vF v model spectrum of Cen X— 4, and the observed Chandra j ACIS-S BI data 
above 0.5 keV. The solid line is the best-fit unabsorbed H-atmosphere plus power-law model 
spectrum with A^h,22 =0.055 held fixed (that is, the intrinsic X-ray spectrum of Cen X— 4, 
prior to absorption by the inter-stellar medium; see Table |l|). The dashed line is the model 
power-law component, and the dashed-dotted line is the H-atmosphere component. The 
two spectral components are equal near ^2 keV, above which the power-law component 
dominates, and below which the H atmosphere component dominates. The H atmosphere 
spectral parameters are r =13. 31^9 km, and kT =0.07 At o^H keV. The power-law component 
is rising in vF v . The crosses are the observed Chandra data, with error-bars in countrate. 
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Table 1. Chandra Observation 
Spectral Parameters (0.5-10 keV) 



H Atmosphere + Power-Law 



a 

-Fx, PL 
kT (keV) 
r 

Total Model Flux 
xS/dof (prob) 

H Atmosphere + Power-Law (Nn fixed) 



a 

-Fx, PL 
kT (keV) 
r 

Total Model Flux 
xl/dof (prob) 



Note. — X-ray fluxes are un-absorbed, in units 
of 10 -13 erg cm -2 s -1 (0.5-10 keV). Assumed source 
distance d=1.2 kpc. 



-, + 0.10 
-0.07 

1.0±0.4 



0.HT 



5 

0.067±0.019 

19 +45 

13 

0.94/9 (0.49) 



(0.055) 

1 o +0 - 6 

1 - u -0.4 

4.6 

074+ 012 
u.u(^_ 005 

13.313 9 km 

11.3 

0.83/10 (0.50) 



Table 2. Cen X— 4 Quiescent Observations 



Ref. 


Obs Date 


kT BB 


X-ray Luminosity 


a 






Instrument 




(passband) 







1980 Jul 28 

Einstein /IPC 
1986 Feb 21 

EXOSAT/ CM A 

1994 Feb 27-28 

ascm/sis+gis 

1995 Aug 16-26 

ROSAT/BRl 
1997 Feb 4-5 

ASCA/G1S 

1999 Feb 9 
SAX/LECS+MECS 

2000 June 23 
Chandra /ACIS-S 



0.32Z 



+0.03 



0.16 



0.13±0.02 



176+ 012 

u - i,D -0.015 



2-3 xlO 32 
(0.5-4.5 keV) 

4-11 xlO 32 
(0.5-4.5 keV) 
4+^xlO 32 
(0.5-10 keV) 

7xl0 31 
(0.1-2.4 keV) 



j+3 



x 10 



1.9±0.3 



2.5±0.5 



(0.5-10 keV) 

(see discussion of § 2.3.1) 
(0.5-10 keV) 

1.7X10 32 \:2 
(0.5-10 keV) 



+0.4 

0.5 



0.066 
0.066 
<0.2 

0.3±0.1 

<0.06 



Note. — The quoted luminosity is the luminosity at the source (i.e. - unabsorbed, using the column density noted), 
source distance d=1.2 kpc; luminosities quoted from references which assumed different distances have been changed. 



Assumed 



References. — 1, 



Van Paradijs et al. 1987 2 



Asai et al. 1996t; 3 



Campana et al. 1997; 4, Asai et al. 1998 



5.E 



ampana 



2t al. 200C ; 6, present work. 



Table 3. H Atmosphere+Power Law Spectral Parameters of Chandra/ ACIS and 

ASCA/S1S+GIS 



(10 22 cm" 2 ) 


a 




kT 
(keV) 


r 

( km ) (rjfes) 


Variable 
Parameter 


Chandra 
Value 


ASCA 
Value 


Xl/dof 


prob. 


(0.055) 
(0.055) 


1.0±0.4 

1 7 +0A 


5 7 +i.o 
°-'-2.4 


076+ - 005 


19+3.0 

1Z -2.4 
12+ 5 


kT 
^X, PL 


0.077±0.007 

r+1.8 

°'°-0.9 


0.085iro°o? 


1.5/47 
1.5/47 


0.01 
0.01 



Note. — 6 Power-law fluxes are unabsorbed (0.5-10.0 keV), in units of 10 13 ergem 2 s x ; Values in parenthesis are held fixed. 



Table 4. H Atmosphere+Power Law Spectral Parameters of Chandra/ ACIS, 
5eppo&4X/LECS+MECS and 4SC4/SIS+GIS 







Chandra + BeppoSAX 












(10 22 cm- 2 ) 


kT r 
( keV ) ( km ) (ofc) 


a 








X^/dof 


prob. 


(0.055) 


0.076±0.007 12.9±2.6 


1.0±0.4 


4.4l^; 9 3 | 


1.9±0.3 


8.6±1.6 


1.0/52 


0.43 



Note. - 



6 Power-law fluxes arc unabsorbed (0.5-10.0 keV), in units of 10 13 erg cm 2 s 1 ; Values in parenthesis are held fixed. 



